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摘 要 : 本 文采 用 分 子 动力 学 方法 模拟 了 温度 363.15 K 和 413.15 KK， 压力 3 MPa 和 6 MPa 条 件 下 ， 硫 化 氧 在 不 同 浓 度 无 机 盐 
溶液 中 的 溶解 度 并 对 Hofimeister 效应 中 反常 现象 进行 了 讨论 。 研 究 表明 ， 相 同 温度 和 离子 浓度 下 ， 压 力 越 高 硫化 氢 的 溶解 度 越 
大 ; 相同 压力 和 离子 浓度 下 ， 温 度 越 高 硫化 氧 的 溶解 度 越 小 ， 相 同 温 度 和 压力 下 ， 离 子 浓度 越 高 ， 硫 化 氧 的 溶解 度 越 小 ， 与 相 
关 实 验 结果 一 致 。 溶解 度 模拟 数据 和 geochem 平台 提供 的 数据 相 比 ,最 小 误差 为 5.5%， 最 大 误差 为 17.7%, 平均 误差 为 13.1%。 
同时 ， 还 发 现 Hofimeister 效应 中 的 反常 现象 是 由 阳离子 的 尺寸 效应 引起 的 。 
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Molecular Dynamics Simulations of the Solubility 


of H,S in Aqueous Solutions 
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Abstract: The solubility of hydrogen sulfide in aqueous solutions is studied by molecular dynamics simulations at 363.15 K and 
413.15 K at pressures 3 MPa and 6 MPa. The notable exception of Hofmeister effect is explained. Results show that the higher pressure 
leads to the higher solubility, the higher temperature leads to the lower solubility and the decrease of the solubility of hydrogen sulfide 
caused by the presence of salt is also observed. The simulation results are compared with the geochem-model data and the minimal error 
is $5.5%, the maximum error js 17.7%, the mean error is 13.1%. It is found that the notable exception of Hofmeister effect is related to the 
cation size. 
Key words: hydrogen sulfide; ion; solubility; molecular dynamics simulation 
人 半 和 工程 应 用 价值 。 
国内 外 的 研究 者 们 做 了 大 量 的 HzS 溶解度 研 究 
H2S 是 自然 界 中 最 常见 的 气体 之 一 ， 广 泛 发 现 工作 ， 得 出 上 百 套 实验 数据 、 几 千 个 数据 点 并 建立 
于 天 然 气 ( 了 J 和 流体 包 里 体 中 41。 随 着 石油 天 然 气 了 理论 模型 外 推 到 实验 所 覆盖 的 温度 、 压 力 范围 之 
勘探 开发 的 不 断 深入 ， 油 气 开采 环境 越 来 越 恶劣 ， 外 ， 但 是 这 些 模型 都 有 一 定 的 局 限 性 。 由 于 Hs 的 
a 剧 毒性 、 腐 蚀 性 以 及 实验 的 温度 、 压 力 条 件 和 成 分 
蚀 ， 溶 解 在 水 基 钻 井 液 中 返回 工地 ， 灾 难 性 的 突 发 的 限制 ， 因 此 采用 分 子 动力 学 模拟 是 一 种 很 好 的 方 
Ce 法 。H2S 在 水 溶液 中 溶解 度 的 机 理 研究 工作 非常 少 ， 
一 些 H2S 含量 较 高 的 气井 站 ,在 我 国 也 相继 开发 了 高 Vorholz 等 ("采用 分 子 模拟 (MC) 的 方法 定性 分 析 了 
含 H2S 的 酸性 天 然 气 藏 , 四 川 盆 地 2/3 气田 含 H2S。 H2S 在 纯 水 和 NaCl 水 溶液 中 的 盐 析 效 应 , 仅 做 了 部 
因此 ,研究 HS 在 水 溶液 中 的 溶解 度 有 重要 的 理论 分 定量 分 析 。Lopez 等 中 采用 分 子 动力 学 (MD) 研 究 
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了 H2S 在 纯 水 中 的 溶解 度 ， 他 们 只 统计 了 气相 区 和 
液 相 区 的 密度 。H2S 在 水 溶液 中 的 Hofmeister 效应 
还 没有 文章 报道 。 

分 子 动力 学 方法 能 够 反映 分 子 间 的 相互 作用 
是 研究 微观 现象 的 有 效 手段 。 本 文 从 分 子 水 平 上 研 
究 HS 在 水 溶液 中 的 溶解 度 ， 并 解释 Hofmeister 效 
应 中 的 反常 现象 。 


1 模拟 细节 


图 1 模拟 盒子 示意 图 


Fig. 1 Simulation cell 


采用 3000 个 水 分 子 和 2000 个 HS 分 子 ， 模 拟 
盒子 的 大 小 由 体系 的 密度 确定 ， 为 了 分 析 问 题 的 方 
便 , 在 分 析 数 据 的 时 候 仅 取 气 相 部 分 LxLxL 的 体积 
来 讨论 ， 如 图 一 所 示 。 模 拟 过 程 中 在 X, Y 和 Z 三 
个 方向 上 均 采 用 周期 性 边界 条 件 。 分 别 各 自 取 108， 
216，324 个 Na ，CL 放 在 模拟 盒 中 ， 配 成 1.94 
mol/kg，3.88 mol/kg，5.83 mol/kg NaCl 溶液 进行 分 
et 水 分 子 模型 选用 键 长 、 键 角 可 变 的 
SPC 模型 ? 。 HoS 模型 用 Shyamal 模型 “， 离 子 模 
型 用 SD 模型 1， 分 子 间 的 相互 作用 力 分 为 短程 范 
德 华 力 与 长 程 库伦 力 ， 分 子 间 相互 作用 势 函 数 为 : 


> 


7CE0 i jl Dy 


ET 人 
G7)=48,[(—) -( + 
ry 广 


其 中 ,oy 和 &; 分 别 为 原子 间 的 尺度 参数 和 势 阱 深度 。 
4 和 gj 分 别 为 原子 i 与 j 的 电荷 , rj 为 原子 i 与 j 间 
的 距离 。&, 为 真空 介 电 常数 。 表 1 给 出 势 函 数 中 的 
各 原子 参数 。 

表 1 势 函数 中 各 原子 参数 


Tablel Parameter of atoms for potential function 


Sigma/ Epsilon/ Charge/ 
Atom 
nm kJ .mol" e 

O 0.3166 0.6502 -0.82 
Hao 0 0 0.41 

S 0.372 2.0789 -0.248 
Hos 0.098 0.0324 0.124 
Na 0.2350 0.5440 1.0 


CT 0.4400 0.4184 -1.0 
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水 分 子 内 的 相互 作用 为 键 和 键 角 围绕 其 原点 的 
简 谐 振动 站]; 


OS h(n) + > so(0-0) O) 
Bonds Angles 
式 中 ,rp 和 06 分 别 为 0.1 nm 和 109.47°*， 甩 和 分 
别 为 465360 kJ/amolnm2 和 386.4 kJ/(mol.rad”)。 
分 子 间 耦 合 尺度 参数 as 和 能 量 参数 & 可 以 根 
据 Lorentz-Berthelot 混合 法 则 [计算 : 


El =A\5E OO =(0,+0.)/2.0 (3) 


为 减少 计算 时 间 ， 势 能 截断 半径 取 0.9nm， 超 
hd 长 程 间 
电力 的 处 理 采用 Ewald 加 和 法 M1。 采用 Velocity 
Verlet 算法 求解 粒子 的 动力 学 方程 ， 系 统 温度 采用 
Nose-Hoover 方法 控制 , 模拟 300,000 步 系 统 达 到 平 
衡 ， 模 拟 统计 参数 300,000 步 ， 时 间 步 长 为 0.8 全。 


结果 与 讨论 


2.1 HS 在 NaCl 水 溶液 中 的 溶解 度 

图 2 给 出 了 363.15 K, 3 MPa, NaCl 浓度 为 1.94 
mol/kg 条 件 下 ， 通 过 模拟 得 到 的 气 液 相 平衡 时 H2S 
和 水 溶液 的 密度 分 布 的 平均 值 。 
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图 2 HS 和 溶液 的 密度 分 布 
Fig. 2 Density profile of H;S and solution 


从 图 2 中 可 以 看 到 气 液 界 面 区 附近 HzS 浓度 比 


气相 主体 浓度 高 很 多 ， 这 是 目前 的 实验 不 能 得 到 的 
数据 。 

表 2 给 出 的 是 不 同 温度 、 压 力 下 H2S 在 NaCl 
水 溶液 中 溶解 度 的 模拟 值 与 geochem 平 台中 查询 值 
的 对 比 ， 最 小 误差 为 5.5%， 最 大 误差 为 17.7%, 平 
均 误差 为 13.1%。 从 图 中 可 以 看 出 ， 模 拟 值 与 查询 
值 曲线 趋势 大 致 相同 证 明 本 文 模拟 程序 参数 正确 ， 
并 且 模 拟 值 都 比 实验 值 要 大 ， 这 可 以 通过 加 入 校正 
系数 来 减 小 误差 。 


表 2 不 同 条 件 下 HS 溶解 度 的 模拟 值 与 查询 值 对 比 站 一 一 一 一 一 一 一 一 一 一 一 一 


2.4 ] 
Table 2 Solubility of H;S in aqueous solutions of NaCl 人 63.15K 
序 误 售 2.0 一 "一 413.1SK | 
T/ P/ 1Nacl/ Ms / mis.exp / 兰 
号 及 差 有 
K MPa 三 16 J 
# mol:kg” mol:kg" mol:kg” % 二 
1 363.15 3 0 1.297 1.102 17.7 1.2 ] 
2 363.15 6 0 2.284 2.03 12.5 
0.8E , ， 1 ， 1 1 1 
3 “363.15 3 1.94 0.933 0.799 16.8 0 1 2 3 4 s BE 
mo (mol/kg) 
4 363.15 6 1.94 1.692 1.469 15.2 ey 
5 363.15 3 3.88 0.737 0.628 17.4 图 4 不 同 条 件 下 Hs 溶解 度 对 比 
6 363.15 6 3.88 1.216 1.153 5.5 Fig.4 Solubility of H;S in aqueous solutions of NaCl 
.1 .8 0.611 0.536 14.0 a 
0 3 3 0 2. 2 Hofmeister 效应 
8 363.15 6 5.83 1.131 0.981 15.3 py 2 ys a 
Hofmeister 按照 不 同 的 盐 类 对 和 蛋白 质 的 溶解 度 
9 413.15 3 0 0.902 0.771 17.0 A L EE 
的 影响 编制 了 Hofmeister 序列 , 离子 的 排序 按照 它 
10 413.15 6 0 1.814 1.629 11.4 ， we Wt eS 
< 一 们 对 溶解 度 从 大 到 小 排列 中 。 通 过 分 析 这 个 规律 可 
es ll 413.15 3 1.94 0.615 0.558 10.2 RS i a 
以 发 现 电 蓓 越 大 或 离子 尺寸 越 小 则 盐 析 效应 越 明 
12 413.15 6 1.94 1.304 1.173 11.2 i ee | 
显 , 显著 的 例外 是 Li 尺寸 比 Na 小 ,但 是 它 的 盐 析 
13 413.15 3 3.88 0.488 0.439 11.2 总 Se A 
效应 却 没 有 Na 明显 。 本 文 尝试 以 H2S 在 无 机 盐 中 
14 413.15 6 3.88 1.002 0.915 9.5 i 本 wp LEN 
的 溶解 度 的 分 子 动力 学 模拟 来 解释 这 一 现象 。 
15 413.15 3 5.83 0.432 0.373 15.8 S 2 Ee E 
为 了 研究 Hofmeister 效应 , 本 文 取 363.15 K, 6 
16 413.15 6 5.83 0.843 0.775 8.7 


MPa, 1.94 mol/kg 时 的 盐 浓 度 来 进行 分 析 。 减 小 Na- 
图 3 给 出 了 413.15K, 3 MPa 条 件 下 , 不 同 浓度 。 ”和 CT 尺寸 参数 15% 和 30%， 其 他 参数 不 变 。 所 得 
的 NaCl 水 溶液 中 HS 的 溶解 度 模 拟 结 果 与 geochem 数据 如 表 3 所 示 : 


I A I Es qe > = a ee 
| 台 查 询 平 台 结果 的 对 比 。 从 图 3 可 以 明显 看 出 ， 表 3 ES 在 NaCl 水 溶液 中 的 溶解 度 
y 一 定 的 条 次 源 的 》 ; 
wy 在 温度 、 压 力 一 定 的 条 件 下 ,NaCl 水 溶液 的 浓度 越 Table3 Solubility of H;S in aqueous solutions of NaCl 
高 ，H2S 的 溶解 度 越 小 。 
EE TT | Sigma+/ Sigma-/ 11INacl/ Mis / 
1.0 上 -| solution E 
4 413.1SK 3MPa nm nm 
Wm 一 一 Simulation Imol . kg mol: kg 
全 一 全 一 h 
Ed 本 NaCl 0.235 -0.440 1.94 1.469 
= 
& NaCl(15%) 0.235 0.374 1.94 1.352 
| ] Na(15%)Cl 0.1998 0.440 1.94 1.205 
NaCl(30%) 0.235 0.308 1.94 1.211 
.4 上 ] 
， Na(30%)C1 0.1645 0.440 1.94 1.361 
0 1 人 
由 /molkg) Na(30%)C1(30%) 0.1645 0.308 1.94 1.187 
图 3 H2S 溶解 度 的 模拟 值 与 geochem 平台 查询 值 对 比 1.sr ; 
a) 
Fig. 3 Solubility of H;S in aqueous solutions of NaCl EE 
er 
兰 1.4 上 
图 4 为 相同 压力 〈6 MPa)， 不 同 温度 条 件 下 ， 日 
不 同 浓度 的 NaCl 水 溶液 中 HS 溶解 度 对 比 。 由 图 | 
可 知 , 在 相同 的 压力 、NaCl 浓度 条 件 下 , 温度 越 高 ， 
Hs 溶解 度 越 小 ， 这 些 与 实验 结果 趋势 一 致 。 
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5 (a) 改变 阴离子 尺寸 大 小 对 HS 溶解 度 的 影响 ; 
(b) 改变 阳离子 尺寸 大 小 对 Hos 溶解 度 的 影响 
Fig.5 (a) Effect of anion size on the solubility of H2S; 


现 


(b) Effect of cation size on the solubility of H2S 

从 图 5 可 以 看 出 减 小 CT 尺寸 , 在 相同 的 条 件 下 
会 减 小 HS 在 盐 溶 液 中 的 溶解 度 ， 减 小 Na 尺寸 ， 
刚 开 始 会 减 小 H2S 在 盐 溶液 中 的 溶解 度 ， 但 是 进 一 
步 减 小 Na 尺寸 ，H2S 在 盐 溶 液 中 的 溶解 度 反 而 增 
加 。 这 说 明 在 Hofimeister 效应 中 的 反常 现象 是 由 于 
Li 尺寸 大 小 引起 的 ， 而 不 是 它 的 自然 属性 导致 的 。 
从 表 3 可 以 看 出 ， 同 时 最 大 程度 地 减 小 Na 和 CT 
的 尺寸 会 最 大 限度 地 减少 H2S 在 盐 溶液 中 的 溶解 
度 。 


3 结 论 


本 文采 用 分 子 动力 学 方法 模拟 了 在 363.15K 和 
413.15 KK， 压力 为 3 MPa 和 6 MPa 条 件 下 ， 硫 化 氧 
在 不 同 浓度 无 机 盐 溶液 中 的 溶解 度 。 通 过 调整 离子 
尺寸 来 改变 离子 大 小 解释 Hofmeister 效应 中 的 反常 
现象 。 研 究 发 现 ， 相 同 温度 和 离子 浓度 下 ， 压 力 越 
大 硫化 所 的 溶解 度 越 大 ， 相 同 压 力 和 离子 浓 / es 
温度 越 高 硫化 所 的 溶解 度 越 小 ; 相同 温度 和 压力 下 ， 
离子 浓度 越 高 ， 硫 化 氢 的 溶解 度 越 低 ， 与 相关 实验 
结果 一 致 。 溶 解 度 模 拟 数据 和 geochem 平台 查询 结 
果 相 比 ， 最 小 误差 为 5.5%， 最 大 误差 为 17.7%， 平 
均 误 差 为 13.1%。 同 时 ， 还 发 现 Hofmeister 中 的 反 
常 现象 是 因为 随 着 阳离子 尺寸 的 减 小 ， 对 盐 析 效应 
的 影响 先 增 大 ， 当 尺寸 小 到 一 定 程度 时 ， 对 盐 析 效 
应 的 影响 反而 减 小 ， 而 这 与 离子 的 尺寸 相关 ， 不 是 
由 离子 的 自然 属性 引起 的 。 同 时 最 大 程度 的 减 小 
Na ，CL 尺 寸 会 最 大 限度 的 减少 硫化 所 在 盐 溶液 中 
的 溶解 度 。 
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